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Introduction. 

§ 1. As a preliminary to attacking the problem oi determining the equilibrimii 
configurations of a rotating mass of liquid, I was led to consider whether some 
method could not be devised for calculating the potential of a homogeneous mass 
in a manner more simple than that usually adopted. What was obviously required 
was a calculus enabling us to write down the potential of such a mass by an 
algebraical transformation of the equation of its boundary, instead of by an 
integration extending throughout its volume. 

There was found to be no difficulty in reducing the calculation to a problem of 
algebraical transformation, but in three-dimensional problems the transformations 
required were, in general, as impracticable as the integrations which they were 
intended to replace. This was because the transformations depended upon a 
continued application of the formula which expresses the products or powers of 
spherical harmonics as the sum of a series of harmonics. 

As soon, however, as we pass to the consideration of two-dimensional problems, 
the spherical harmonics may be replaced by circular functions of a single variable. 
The transformation now becomes manageable, and for this reason the present paper 
deals only with two-dimensional problems. 

The first part of the paper contains a short sketch of a theory of two-dimensional 
potentials. I have, however, confined myself strictly to such problems as are 
required for the solution of the main problem under discussion, namely, that of the 
rotating liquid ; the method does not attempt to be one of general applicability. 

The Potentials of Homogeneous Cylinders. 

General Theory. 

§ 2. We shall suppose the cross-section of the cylinder of which the potential is 
required, to be bounded by a single continuous curve S enclosing the origin. Let 
(322.) K 2 25.IL02 
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V be the potential of this cyhnder, supposed to be composed of homogeneous 
matter of density p. 

The value of V must be finite and continuous at all points except infinity, and 
its first differentia] coefficients must also be finite and continuous at all points. 
Also V must satisfy V"^ V ~ at all points outside S, and V^ V ^ — 47r/3 at all 
points inside S. At infinity V must vanish, except for a term proportional to log r. 
' These conditions suffice to determine V uniquely. For if there were two distinct 
solutions V and Y\ the function (V — V) would satisfy V^(Y — V^)-^0 at all 
points of space, would be finite and continuous, together with its first differential 
coefficients^ at all points of space, and would vanish at infinity, except for a term 
proportional to log r. The only solution satisfying these conditions is known to be 

Y -».. y ::::z: 0, heuce auy function V satisfying the conditions laid clown above must 
be the potential of which we are in search. 

§ 3, Let us use polar co-ordinates r, 6 in conjunction with orthogonal co-ordinates 
X, y^ and let us also introduce complex variables f, 07, defined by 

Let us suppose the eciuation to the curve S ti^ be written in the form 

f{t v) = o (1), 

and let us imagine that this equation is solved explicitly fbi' f in the form 

£-F(^) (2). 

In general 1^ {rj) will be a multiple-valued function of r], and the equation (1) 
may, and probably will, be satisfied for other vidues of f and r] than those which 
occur on the surface S. 

Let us, however, suppose that we have succeeded in finding one value of F (r)) 
such that this value is a single-valued function of t] at every point of S, and is 
equal to f. Let us suppose that we have succeeded in expanding this value of F (7^) 
in a series of ascending and descending powers of r], these series each being supposed 
convergent at every point of S. Let us write 

7 7 

^ iv) == «0 + ^iV + ^^ + . . . (4), i/i (tj) r=: ^ + ^1 + . . . . (5). 

We shall consider only the case in which the surface S has the plane 6^ = as a 
plane of symmetry. Li this case the equation f {^, rj) ^=:=: remains the equation to 
the curve after the sign of 6 is changed, and we therefore have as a second form 
of this equation, 

There is therefore a solution of this equation expressing 7; explicitly as a function 
of ^ in the form 

rj = t {^) =: 4> {^) -r ^|J {^) (7) 

wliere F, </>, \p are defined by equations 3, 4 and 5. 
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Since we have supposed equEition (2) to represent the complete solution of 

equation (I), the equation 

^ = F(a (8) 

must necessarily represent the complete solution of equation (6), if the meaning 
of F remains unaltered. Let us suppose that ilo values of the multiple-valued 
function F (cf ) (or, what is the same thing, F {x + iy) ) are exhibited on the appro- 
priate Riemann's surface. The locus of points at which this function is equal to 
X '— iy (also a definite and unique function of position upon this Riemai^n's surface) 
will be the complete system of points satisfying equation (6), and will therefore 
include the curve S. This curve will not, however, be drawn upon a plane, but 
upon a Riemann's surface. 

Now so long as the curve S does not possess a cusp or branch point we have at 

every point of S, 

d¥ (|) __ d (x - ill) 

d^ '"' d (x '{- iy) 

and hence it follows that no branch point of the Riemann's surface can lie on the 
curve S. Let the curve S be a single closed curve surrounding the origin, and it 
follows that we can always so arrange the Riemann^s surface that no branch line 
shall intersect the curve S. In this case it is possible for a solution of the form of 
equation (6) to satisfy at every point of the curve S, the curve now being regarded 
as a closed curve lying on one sheet of a Riemann's surface. By simply interchanging 
the axes of y and — y we can imagine the function F (t;) represented on the same 
Riemann's surface, and we see that at every point on the old curve S we shall have 

^ = F (^) = <^ (t,) + rMr;). 

We have therefore proved that if a particular solution of the form of equation (3) 
can be found to represent the curve S, then solution (7) will represent the same 
curve. We have also seen that if a family of curves is described by continuous 
deformation starting from some curve Sq, the general solution can be arrived at by 
continuous variation starting from the solution for the curve Sq, so long as no 
member of the family j)ossesses a cusp or branch point. For our present purpose 
it is sufficient to have proved that if we have, at every point of the surface S, 

f = </. (t?) + i// (.,) (9), 

then we have also, at every point, 

.j:-<^(0+r>(a (10). 

Let us now introduce a new function ;)( defined by 

where C is a constant^ 

We have by differentiation, 

^x/^^ =- <t^ iO + i^ iO - V ■ ■ (12), 3x/ai? = <^ W) + tH^) - f • • (18), 

and therefore at the surface 8, by equations (9) and (10), 9x/^f ^^ ^x/^V = ^- 
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It follows that X ^^^ ^ constant value at the surface S, and this, by a suitable 
choice of the constant C, may be taken to be zero. Also it follows that d^jdn = 
at all points on the surface S, where djdn denotes differentiation with respect to the 
normal. We therefore have at every point of the surface S 

X=0 . (14), dx/dn = ...... . (l5). 

§ 4. Let us denote the potential at a point inside the surface S by V,-, and that at 
a point outside S by Yq. Let us examine, as a trial solution for V, 

Y,= ^p\c+\'cl>(i)di+\"<j>{7,)drf~-^yj] (16), 

L •'0 -'0 J 



^ (0 <^'i + \ "P (v) dv \ (1^)- 






Since r) = re~'^, the greatest value of | ^ | at any point on S is equal to the greatest 
radius, say R^, which can be drawn from the origin to S. Since the series (j) (t)) is, 
by hypothesis, convergent for all points on the boundary, it follows that the radius ot 
convergence of the power series (j) (rj) must be greater than H^, and hence that (j) (r)) is 

convergent at all points inside the surface S. Hence also (f) (tj) cIyj must be convergent 

at; all points inside S. The same is obviously true if i is written instead of rj, Hence 
it follows that if V^ is defined by equation (16), then V^ and its first differential 
coefficients will be finite and continuous at all points inside S. 

Li a precisely similar manner it can be shown that if Vq is defined by equation (17), 
then Vq will be finite and continuous at all points ou.tside S except at infinity, and 
that the first differential coefficients of Vq will be finite at all points outside S. 

From what has been said it follows that Vq and V,- are finite at the boundary. 
The value of V^- — Vq at the boundary is wpx^ and this vanishes by equation (14). 
Hence V is finite and continuous at all points of space (except infinity). 

Since the series <^ (r)) and xjj (t)) have been supposed to be convergent on the 
boundary, it follows that the first differential coefficients of V must be convergent on 
the boundary, and hence that these differential coefficients are finite at all points ot 
space. At the boundary, 

dYi/dfl •— ■ dVQ/dn = Trp dx/ofi =0, 

by equation (15). Hence it follows that the first diffei^ential coefficients of V are 
finite and contimtous at all points of space. 
At a point inside S, 

V-^V = V^^V,- = 4.TTP ^ , V j G + f' <^ {^) d^+ f V (^) d^ - ^7]\ ^ - 471/), 
and similarly at a point outside S, 



'^Q r f 00 , CO 



v^V = v^Vo = 477P ^^ I _ ^ ,|; (f ) ./^ + _ ^ t/, {-n) d-n^ = Q. 
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Lastly, by actual integration, we find as the value of Vq from equation (17), 

Vo = vrp I - 6, (log ^ + log ^) + h, Q + I) + i&3 {j. + |-) + . . .} (18). 

The term — h^ (log £ + log rj) may be replaced by — 26^ log r ; hence V vanishes at 
infinity, except for a term proportional to log r. 

We have now seen that V satisfies all the conditions which must be satisfied by a 
potential ; hence the potential willbe given by equations (16) and (17). 

§ 5. Calculated by direct integration, the value of Vq is 



■Vo - - 



log liV drUl9\^ 



where the integration extends throughout the cross-section of the cylinder, and R is 
the distance of the point r, 6 at which the potential is being evaluated, from the 
point r'', 6^ of the cylinder. We have 



E^ = r^ + r'^ - 2i¥ cos {9 - 6') = (r - rV^<^-^')) (r - r'e~'(^--^'^y 



Now if 



> I r^ I , we have 



log (r — r^e '^^ ^'^) = log r + log ( 1 
and hence 






log 7' 



rW^' J I'/e^^Y 



re^^ 



\ Q^P^^ 



f^j^e' r's-^'X 



lo2^ K^ = 2 loff r — — .. A 



■ie 



w J 



• t • 



Upon integration we obtain 



V 










p 2logr j r'dr'dff — ^ — 



r'-^ti'"' d7'' d6' + ^'-ll r'-^e-'O' di-' d6' 

7) J 



(19). 



The cylinder being symmetrical about the plane # = 0, we have 



/* f 



r'-Wdr'dd' 



and hence equation (19) becomes 



V V 



r%~''^' dr^ d6^ 



v^'^ cos ffdr^d&. 



Vq = — p (log f + ^og yj) r^dr^dd^ + p {~^~ + ~~" ' 



r^^ cos 6' dr^d6^ + • 



(20). 



Let the area of the cross-section be A, and its centre of gravity be at x = a, 
so that 



XX — ■ 



rUlr' dS' 



XX 



a 



r^'^ cos 6' dr^ dO ; 



then equation (15) becomes 



/I 1 \ 
Vq = — pA (log i + log 7?) + pAa { y H ') + ••.. 

Comparing this with equation (18), we find that 

A = Trbi ........ (21), Aa = irb^ ........ (22). 

These last two equations enable us to find, by a process of algebraical*! transfor- 
mation, the cx^oss-section and the centre of gravity of the cylinder, whose bbundary is 
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given by equation (1). 1 have made no attempt to discuss these results from the 
point of view of pure mathematics ; we therefore pass at once to some simple 
illustrations of our theory. 

Circ u lav Cy lin der\ 

§ 6, Let us attempt to find the potential, cross-section, and centre of gravity of the 
cylinder^ of which the cross-section is the circle 

(x — 6*)-'^ 4- y^ =^ a°' . 
or^ in polar co-ordinates, 

r^ z=z afi -|~ 2er cos 6 -— a 

The ^5 7] equation to this curve is 

^ri = a^ »j-. c (1^ 4" ^) "^ 
or, solving explicitly for £ 

f — c + c^jiji — c) 



( • > ^\ 






1 ,mim "!!J\ 



c- 



(^2oj, 



Now the minimum value of 



7] 



which can occur on the curve is a — c, and this is 



TOater than c if a > 2c. In this case we have 



Ti — G 

and therefore the solution for 'c is 






+ -~ + 



t — 



a' 



d" 






c 



^M-'^ T'^' 



a? 






V 



1" . .3 ""r 



4 « 



We accordingly have 



— ITp -^ C (4* 4" '^' 



V; 



TTP 



{ 



•1 



1 \ 



(^ log |i? + ca^ ( "V" 4" —- ) + i 



^' 



TTCf^" 



.0. 






1 -^ 






1 \ 



• 6 \ ,«uJ W /• 



rrcci?. 



If we determine C from the condition that Vq and Y- shall become equal at the 
boundary, we have the known values for Vq, V,,. The equations for A and a reduce 
to A = na^, a =: c, and these, again, are obviously in agreement with the 
known results. 



a e <> 9 <} 



Elliptic Cylinder, 

§ 7, Let us next find the potential of the elliptic cylinder 

ax^ 4- 'by'^ = 1 
The f, 7) equation to the surface is 

{a - h) {i- + rf) + 2 {a \- l>) iv '"^ 4, 
or, solving explicitly for (;, 

1 



a — b 



(a 4- &) '^i 4" y ^icihrf 4~ 4 {a .— b) 



{^jU I A 



• • - • - I amJ t/ /• 
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At points on the boundary the mhiimnm vakie of 4a& | -^^ j is 4a ; hence, provided 
a^ > (a — 6)^, a convergent expansion for ^ is 



^ 



a 



a — h 



(« + ?')^ + 2v/(a?>)r?(l + i ^^ 



1 



(« - 1^ 

7 t 



* « • 



1 



a — & 






• ♦ • 



(30). 



From this expansion we obtain at once 

27rp 






^a H- ^/6 



{^ax^ + v^6r) 






a 







, (olh 



TT 

We can obtain Vq in series at once ; if we require its vahie in finite terms w^e may 
proceed as follows : — - 
From equation (30), 

1 , a — h 



-co / 










77/) 



tj V V 

.CO 



% - -y-j; T I + • • • ) ^^'^ + same function of £ 



= 7rp 



>? 



a 



— (a + ^) >? + \/ ^(xhyf + 4 (a — 6) I + --, f --^ ->; 



drj 



4" same function of ^ 



irp 



V 



a — 1) 



^{ah) rf + „ y/4:abyf + 4 (a — 6) 



+ '~77T\ ^^E (2v/(<^^^) 1 + y ^cibrf + 4 (a — 6) 

+ same function of £ . . . (32). 

The results obtained agree with the known results if ^a and ^h are taken of the 

same sign. There is a second solution, obtained by changing the sign of one of these 

roots, and this corresponds to a mathematical ellipse of which one axis is negative. 

The full significance of this will appear later. 



Ewpansion in Potvers of a Parameter, 

§ 8. When the equation to the surface is of a degree higher than the second, it will 
not, in general, be possible to obtain a solution in finite terms of the form of equation (29). 
Suppose, however, that the surface forms one of a family of surfaces, the family being 
described by the variation of a parameter c, and let this family be chosen so that the 
surface c = is one for which the complete solution is known. Then as we know 
the value of f when c = 0, we shall assume that it will be possible to find the general 
value of ^ in a series of ascending powers of c, and the equations determining the 
various coefficients of c will have a unique solution. 



VOL. CC— A. 



n 



74 MR. J. H. JEANS ON THE EQUILIBRIUM 

The solution for ^ will be a function of rj and c. The maximum values for j rj | and 

1 . 

— will therefore be functions of c, and will be expressible as series of ascending powers 

of c. Substituting these values respectively in ^ (rj) and ^ (rj)^ we obtain two series 
in ascending powers of c. Let the radii of convergence of these two series be c^ and c^ 
respectively. Then <^ (tj) is convergent inside and over the boundary provided c < c^ 
and \jj (tj) is convergent outside and over the boundary provided c < c^. The solution 
for f found by expansion in powers of c will satisfy the conditions which it is assumed 
to satisfy in § 3 so long as c is less than the smaller of the two quantities Cj and Cg. 

§ 9. Let us suppose our solution in powers of c to have been obtained. The values 
of V,- and V^ can be regarded as power series of the variables f or rj, and as such will 
have circles of convergence in space, having the origin as centre. Let V,- be 
convergent inside a circle of radius R^, and Vq outside a circle of radius R^. When 
c = 0, R^ < R^ ; the circle R^ being wholly inside, and the circle Rj being wholly 
outside the surface S. When c reaches the value Cg, the circle R^ touches the surface 
S ; and for values of c greater than e^, the circle Rg either intersects, or lies outside, 
the surface S. 

Suppose that Cg < Cj : then for values of c such that 

(V J, 

the circle Hi lies wholly outside S, while the circle R^ does not lie wholly inside S. 
Thus if V^ and Vq are defined by the same power series which have been found for 
them for values of c less than c^, V^ will be convergent at all points inside S, but Vq 
will not be convergent at all points outside S. We have now to inquire whether it 
is possible for Vi to give the true value of the internal potential, even when the series 
found for Vq fails to represent the external potential 
Let the equation to the surface be written in the form 

y (f , ij, c) =:: . » . . . . . . , . (33), 

and let the solution be written in the form 

f = ^(t^, c) + t/i(?7, c) . . . . . . . . . (34). 

We shall only consider the case in which c < c^, so that our former function ^ (t)) is, 
by hypothesis, convergent at the boundary and at all points inside it. We take 
(j) {t], c) to mean the same thing as our former ^ (?;), and suppose \b (t;, c) defined by 
equation (34). Now the value of i cannot become infinite at any point on the 
boundary, since we suppose the bounding surface not to extend to infinity for any 
value of c. Hence it follows that our function \jj {t], c) will be finite at all points of 
the boundary so long as c < c^ Moreover, in the region in which our former \jj (rj) is 
convergent {ie., in the region r > R^), xfi (t;, c) must become identical with \fj (rj), and 
will therefore be finite. The function 

ijj (t), c) + \Ij (£ c) 

will therefore be finite at the surface S, and will vanish at infinity to an order at 
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least equal to 1/r. This fimction is real at every point of space, and is a solution of 

Laplace's equation. Hence the function is finite at every point of space. 

It is therefore clear that if, in equations (16) and (17), we replace ^(r\\ ^{rj) by 

(j) {t}^ c), ^ (t), c), and make the corresponding changes in ^ (i), xff {^), we shall have a 

solution which satisfies all the conditions of the problem, subject to the single 

condition that c is less than c^. At infinity xjj (17, c) will be capable of expansion in 

the form of equation (5), 

xj; (rj, c) = hjj) + hjr)^ + hjrf + . . . , 

and we now see that equations (21) and (22) can be obtained in the same manner as 
before. 

It is therefore clear that the values obtained for V^:, A, and a will be the true 
values, even if they have been obtained by the use of divergent series, provided only 
that the series for V^- remains convergent up to the boundary. 

In the case in which c^ < c^, a similar proposition is true for values of c such that 

§ 10. \A' e now consider the case in which c is greater than either c^ or Cg. Let the 
values of V^ and Yq which have been found for values of c less than either c^ or c^ be 
extended, by a process of '' continuation," to points outside their circles of convergence, 
and let the values so obtained define the functions V^- and Yq. These functions will 
have certain infinities, the position of these infinities depending upon the value of c. 
When c = 0, all the infinities of V,; lie outside S ; all the infinities of Vq lie inside S. 
Let us suppose that up to some value of c, say c = Cg, no infinity crosses S, but that 
(if possible) at c = c^ one of these infinities is found on the boundary. The values of 
V,- and Vq are functions of c, and V^- — Vq satisfies the requisite algebraical equations 
from c = until c is equal to the smaller of the values c = Ci or c^. Hence it must 
continue to satisfy for all values of c, until the condition found in § 3 is violated, ^.e., 
until the value of c is such that the curve possesses a cusp or branch point. Also V^ 
and Vq satisfy the requisite conditions of finiteness, uniqueness, and continuity until 
c reaches the value Cg. Now as c approximates to Cg from the direction in which 
c < Cg, the value of V at same point of the boundary (viz., the point at which the 
infinity occurs when c = Cg) will increase indefinitely, becoming ultimately infinite 
when = Og. This value of V will, however, give the true solution for all values of c 
less than Cg, and there will be a superior finite limit to the value of V. It therefore 
follows that there can be no value Cg at which an infinity crosses the boundary, and 
the values of V^- and Vq found by '^ continuation '' of the power series will give the 
true values of V^ and Vq until the whole solution is invalidated by the occurrence of 
a cusp or branch point. 

Summing up, it appears that we may neglect the question of convergency of series 
altogether : so long as the values obtained for either V^ or Vq are possible values, 
they must be true values. But care must be taken not to pass through values of the 
parameter such that S possesses a cusp or branch point for these values. 
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§ 11. Illustrations of these remarks are afforded by the examples of §§ 6 and 7. In 
ecjnations (23), (24), and (25) let us regard c as a variable parameter, so that, as c 
varies, the equations represent the different members of a family of circular cylinders. 
We have seen that the series obtained for f only remains convergent so long as c < -|-a. 
Equations (26) and (27) represent the values of Vi and Vq expanded in powers of c. 
Now the value of V^ is convergent, no matter how great c may be, hence we know 
that this represents the true value of V,- for all values of c. The value of Vq has as 
its circle of convergence the circle r = c, and this intersects the boundary as soon as 
c attains the value c = |-a. Hence the value obtained for Vq will fail to give a true 
solution as soon as c exceeds the value c = |a, although it will always give the value 
of Vq at points outside the circle r =^ c. Again, in § 7, let us regard ^/h as a variable 
parameter. The value found for V^ is convergent for all values of s^/h, and therefore 
represents the true solution for all values of ^6, provided that we start from a true 
solution, and do not pass through a value of ^h at which a cusp occurs. Under this 
same condition the series for Vq will give the true value of Vq at all points at which 
it is convergent, and the expression given in equation (32) will give the true value at 
all points, this being the expression for Vq which would be found by *' continuation " 
of the series, or (what is the same thing) by the methods of § 9. 

The ellipse does not possess a cusp except in the critical cases of ^h = 0, -^/6 = oo . 
In the former the ellipse reduces to a pair of parallel lines, and the points at infinity 
rank as cusps. In the latter the ellipse reduces to a doubled straight line joining the 
points ^ = lb ^"^^ and, again, these points rank as cusps. Hence a solution will- 
remain the true solution so long as ^/b does not pass through either of the values 
^h z= or CO , i,e,, so long as A^h does not change sign. This is the meaning of the 
condition found in § 7, that ^^a and ^Jb must be taken with the same sign. 

We can see this from another point of view^, as follows. The solution (29) can be 
exhibited on a Riemann's sin*face of two sheets, the branch points being given by 

abif •=. b — a. 

When a = b (?.6., when the ellipse reduces to a circle) these points coincide in the 
origin, and destroy one another. As b increases, the two points move along the axis 
oi X, and ultimately meet the curve when b = oo at the points a? r= -^ a~% i.e., meet 
the curve at its cusps as soon as cusjds occur. Similarly, as b decreases from the value 
b = a, the branch points move along the axis of 3/, and meet the curve when 6 = 
at the points it :^ A- m , 

Deformed Circular OyUnder. 

§ 12. As an example of expansion in a series of powers of a parameter, let us 
consider the cylinder of which the equation is 

T^ := ci? 4" 2cr^' cos nO ......,, , (35), 

or, in ^, 7) co-ordinates, 

^77 =^ a^ + c {^'^ + 7]'') . (36). 
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We are in search of a solution for f expanded in powers of c in the form 

£ ■=. u^ 4- u^c + u,^c'^ + ii^^c^ + 

in which u^^ u^^ . , . are functions of t]. 

Substitute this assumed solution in equation (36), and we obtain 

=: a^ + c {tj^^ + tf/ + n^o"-! u^c + '}pi {n — 1) t^,/^"-^ t^^^ c^ 
+ nu^''^^ 11^ c^ + -|-7^ (71 — 1) [n •— 2) -z^q''""^ t^j^^ c^ 

+ 71 (n — l) 7i(^''~^ t^^ t/^c" + UIIq'^^ Uq c^} 

+ terms of degree 4 and higher in c. 

Equating the coefficients of the various powers of c, we obtain 

rjK^ — a^ 7]i{^ = 7)^' + Vq% 

rjU.2 =^ nn^''"'^ ifp 7]iL^ = \n {n — 1) u^'''^u{' + nii^'"'^ ii.^^ 
7)11^ = ^n {n — 1) {n — 2) i^q^'"^ t^^^ + 7^, (n — 1) iiQ^'^^itjih + ^^^^^'"^^^3? &c. 
Solving these equations in succession, we obtain 



u 










w 



V + 



a- 



In 



V 



?2 + l' 



ic^ z=z na 



2u-^2 



1 



+ 



./j+i /' 



t{^ 



\n{n-' l)a'''^^r)''^^+ - 



n {2n - 1) 



V 



n+l 



a 



4«^4 






7] 7}' 

01 (3?i — 1) a^"'~* 



t^' 



4< 



= 1-n (71 -- 1) (71 - 2) a'^^^-S'""^ + ^I'J'^^J.^ c&~^ 



9. 



V 



+ terms of lower degree in r/, &c. 



Hence we obtain at once 
V, = C + rrpl — £^ + 



^n _|_ ^u 



n 



c + i (n - 1) a?^-^^ {^^ + ^^^) c^ 



(37), 



+ iV (^^ - 1) (^^ - 2) a^''~^ {e'" + V^") ^^' + . • . 1 
=: C + 7r/)j - 7- + J~ + (ri — l)a^^'~H^] ?'" cos 7i^ 

+ (^^ •- 1) (^^ -^ 2) cr'^-^ c^^^^^ cos 2ne + ...}. . (38), 
A = TT {a'- + nar"-^ c^ + f ?i (n ««. 1) a*''^^^ ^t + . . . }^ 

and, except in the special case of 7i = 1, a = 0. 

In this way we can, when c is small, write down the potential to any required 
degree of accuracy. 

Deformed Elliptic Cylinder. 

§ 13. As a final illustration, we shall find the potential produced by a small 
deformation of the surface of the elliptic cylinder 

^r^ = a^ + a^{e + rf) (39). 

Let the deformed surface be 
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CO 



#/7 = a^ -I- a^ (f 3 _j- ^^^) + 2 &„ (f " + 7?") . ..... (40), 



and let the solution be, as fai- as first powers of &'s, 



CO 



t ^ p \^ ]. p /4I \ 

1 

where £, is a function of 77, and ^q is the sohition when all the Vs vanish. 

Substitute solution (41) in equation (40), and equate coefficients of h^, and we obtain 

Lr) = 2c% £o£, + ^^' + 7]'\ or, solving for ^„ 

L = do" + i') I (v — 2% £o)- •"«••' • (42). 
We can express ^q in the form 

^Q = a?; •+ ^/r] +••'•••••• (43)» 

Hence equation (42) becomes 



f; 



b i^i 



r) (1 — 2a2a) — 2^3 ^77 ^ + . . . 

.n / 



1 •— 2ao5t ' \l~-2ao«. (1 — 2aoa)^ / ' 



Hence we find as the value of Yf, 
V,- = 7Tp {C »^ £7? + |a (^2 +^s)} 



11 = 00 



^,,^1^ '^ [1 — 2cL2a n V 1 " 2%a {1 — 2ac^af / 7'^ 

1 

4" terms in 6,,+4,5 6,,^g, &c. . . « j> , ., . . (44). 

The value of ^q given by equation (43) is a solution of equation (39). Sub- 
stituting this value, and equating the coefficients of the two highest powers of 77, we 
find as the ec|uations determining a and ^, 

a ~ ^2 (1+ a^), ^ =: a^ + c% {2a^ + 1), 

equations which will be required later» 



EoTATiNG Liquid Cylindeb,. 

Geneoxd Theory. 

§ 14. We now pass to the main problem before us, and consider the equilibrium 
of a cylinder rotating with angular velocity m. 

The equation to the cylinder for a rotation equal to zero is 

When the rotation co is different from zero, we shall suppose the equation to the 
surface referred to its axis of rotation as origin to become 
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iv ^ '^^ + «i {^+ v) + ^3 (1^ + ^^^) + • • (46), 

or, in polar co-ordinates, 

r^ — a^ + 2a^r cos 9 + 2%^^ cos 29 -\- . , (4^). 

This equation is not sufficiently general to represent all cylinders v/liich are 
symmetrical about the initial line. The value of r"^ being known at the boundary, 
we shall always be able to find a function v such that v is finite and continuous, 
together with its first diff*erential coefficients, at all points inside the boundary, 
and such that V^v = inside the boundary, and v = r^ at the boundary. The 
value of V near the origin can be expanded in the form 

I) z=z a^ -|- 2a3^?' cos 9 + 2a^r' cos 29 -{- (48), 

and this series will have a circle of convergence, say r = R. It is only when the 
curve lies wholly inside this circle that the cylinder can be represented by an 
equation of the form of (47). 

Let us, however, attach a conventional meaning to equation (47) in the case in 
which the right hand becomes divergent at the boundary, as follows. I.-et us suppose 
that the value of the function v given near the origin by equation (48) is calculated 
from its known values inside the circle r = R, the values outside this circle being 
obtained by a process of " continuation.'' Then we shall suppose equation (47) to 
represent the locus of points at which r^ =:: v. 

Obviously, with this convention, equation (47) is sufiiciently general to represent 

any surface. If this surface is to give an equilibrium configuration under a rotation 

0), we must have 

V-i +;! 0)%'" ::::= a constant (49) 

at the surface. Now V,- + Trpr^ is a spherical harmonic at all points inside the 
surface, and equation (49) can be written in the form 



0)2 \ 



(Yi + TTpf'^) — • TT/) ( 1 — - - — j r^ == a constant, 



0)3 \ 



or, what is the same thing, 

(V^ + TTpr"^) — TTpil — ^r" j V z=z: a coustaut .' , . . . (50). 

This equation is satisfied at the surface S, and each term is a solution of Laplace's 
equation at every point inside S ; hence the equation must be satisfied at every 
point inside S. 

Now V,- can be calculated by the methods already explained, and we obtain an 
equation of the form 



n — \ 



which gives the value of V,- at all points inside a certain circle of convergence. The 
value of V inside its circle of convergence is, from equation (48), 



V = a' -f t a„ {t + yf) ; 
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hence equation (50), at points inside both, ciixles^ becomes 

— Trpcfr" 1 1 - 



0) 



9 



7rp 



')l= 00 



■^ rrp t \fu{<^^i, (h^ , , ,) 



a 



II 



1 



n^\ 



9 



7r/3 / 



(^" + '>?'') — a constant, . (51). 



Ct) 



27rp 



(52). 



We must therefore have, for all positiYc integral values of n 

The condition that equation (46) or (47) may represent an equihbrium configuration 
is therefore 

and the value of cu" is sucii ttiat each fraction is equal to 1 



27Tp 



. . » ( o L 

Since for any 



equilibrium configuration the axis of rotation must coincide with the centre of gravity, 
we see that if the coeflicients of (46) satisfy (53), the curve must be referred to its 
centre of gravity as origin. 

The points of bifurcation are given by the Hessian of this system of equations. 
For our purpose this may be most conveniently written in the form 






1 - f ) 

2lTpj 















1 



or 



\ 



o 



LTTp! 



. (54). 



Our method will, for reasons already explained ^ break down as soon as a cusj) or 
branch point occurs on any linear series. Any solution will be a true solution 
provided we can pass from this solution to another, known to be a true solution, over 
a path through a system of linear series, without passing through a point at which 
a cusp or branch point occurs. Now the occurrence of a cusp or branch point indicates, 
in the physical problem, the division of the mass of fluid into two separate masses, 
and when this occurs the solution breaks clown for a second reason also : for 
equation (49) is only true when the surface of the fluid is continuous ; when the 
surface consists of two distinct parts, the constant on the right-hand side has different 
values for the two parts of the surface. 

This limitation will not cause trouble in the present investigration. For we are 
only desirous of tracing the changes in the confi-guration of the fluid up to the 
separation into two parts, and even if our method had enabled us to proceed beyond 
this point, it would have been fruitless to do so. 
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The Series of Circular Cylinders. 

§ 15. If it were possible to calculate the /'s and solve equations (53) in the most 
general case, we should arrive at a complete knowledge of the system of linear series 
of equilibrium configurations. This being impossible, we shall start from a known 
series, and calculate successive series by determining the various points of bifurcation. 

Now we know (§ 12) that for the small values of the a's, f^ is of the form 

fni^i, ^0, . . .) = -^' . Hence there is a solution of the system of equations (52) 

given by a^ =: c/g ==...== 0. 

This is the series of circular configurations, and corresponds to the series of 
Maclaurin spheroids in the three-dimensional problem. When co^ > 27rp the solution 
breaks down physically, since the pressure at every point of the liquid becomes 
negative. In fact, when m'^ reaches the value o)^ = 27rp the series gives place to 
a series of annular forms, for each of which o)^ has the critical value. We can adjust 
the radius of the annulus so as to give any desired amount of angular momentum 
pxeater than the critical value which occurs in the circular confip^uration when 

CO^ zzz 27Tp, 

Points of Bifurcation on Circular Series. 

§ 16. To searcli for points of bifurcation on this linear series, we replace f by a^/n 
in equation (54). Every term in the determinant on the left hand now vanishes, 
except the terms of the leading diagonal, and the equation reduces to 



'> \ 



1 - ( 1 - ■'"■ 



< 



» \ 



2 ""(^ "^sttJ; • • • 



1 / 0)2 \ _____ 

n \ 2'Trpj 



The different roots correspond to the different integral values of n, and are 

given by 

n = 1, 2, 3, 4, 5, ... CO, 



0) 



2ii 



'/2ttp = 0, -5, -666, 75, '8, . . . 1. 



The first point of bifurcation (n = 1) may be rejected at once, the critical 
'' vibration'' being merely a displacement of the entire cylinder as a rigid body. 

A displacement in which a,i, only occurs for the single value n -=. s will alter the 
potential energy only by a term proportional to the square of cv Hence the 
principal vibrations correspond to the different values of n from 2 to oo , and are such 
that a« only occurs for a single value of n in each. When w^ = 0, all these vibrations 
are stable. When oy^ reaches the point of bifurcation of order 5, the vibration of order s 
becomes unstable, and, since there is only one point of bifurcation of order s, this 
vibration remains unstable for a.l values of o)^ greater than the value at this point of 
bifurcation. We therefore see that by the time that (J' reaches the limiting value 
271/), every vibration is unsiable. 

VOL. cc. — A. M 
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lie Series of Elliptic Cylinders. 

S 17c The stable linear series is that of order 2, Tlie initial deformation is therefore 
elliptical, and the point of bifurcation occurs for the value cd^ ==; irp. 
The ecjuation to the surfiice is initially 

^ =r. a^ + a^ (P + 7|S) . , , . . . . . (55). 
In § 13 we found (equation (44)) as tiie corresponding value of Vi 

y)^7rp{G~~i7^+ia{e+7f)} ■ (56), 

in which a is a root of 

and this is true however great a.^ may be. The values of /J^. (O, a^, 0, 0, . « .) can 
accordingly be written down at once. We have 

and all the other functions vanish. 

Hence there is a general solution of equations (52) in which, all the €^;/s vanish 

except a.,, and 

J- '^ . 

|-a ==: (l — (xyf^TTp) Cr'o « , . » ^ . . . (58), 

where a is given by equation (57)« 

This is the linear series of which we are in search. It is obviously a series of 
elliptic cylinders, and corresponds to the series of Jacobian ellipsoids in the three- 
dimensional problem. From ecjuations (57) and (58) we have 

(xP-'/iTp = 1 — a^ , . . . . . . . » (59), 

We therefore see that as we move along this series the value of or continually 
decreases from irp to 0. The angular momentum, however, increases from a finite 
to an infinite value. 

The Remaining ( Unstable) Series, 

§ 18. Before searching for points of bifurcation on this series, let us briefly examine 
the series passing the other points of bifurcation on the circular series, these series 
being known to be all unstable. Near the point of bifurcation the form of tlie series 
of order n is 

^r) =: a^ '~\- a,, [t' + ri'') . . . . » . . = (60). 

In § 11 we have calculated the values of j^. (0, 0, . . . a,,, 0, , . .) as far as a,/. If 
we neglect a/, it appears that all these functions vanish except /„ and that /,, is 

given by 

/; (0, 0, . . . a,,, 0, . . .) = a^jn + I {n — 1) a?''~"'a,?. 

The series is accordingly given by equation (60), until a/ become appreciable. 
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The value of Ci)'^ is given by 

1 — myZrrp — I In + 1 {n ~~ 1) (J^^''-%,^ 
It therefore appears that the angular velocity decreases as we recede from the point 
of bifrircation. 

Throughout the series the bounding curve is periodic in 9 with a period 27r/;?, and 
it is easily seen that at the remote end of the series is a curve consisting of n equal 
and symmetrically arranged arms, these arms extending in the limit to infinity, and 
being of zero breadth at all points except in the immediate neighbourhood of the 
origin. 

In. fig. 1 I have drawn the cin^ve for the case of n = 3, a^ 



1 



The 



error m 



the value of V,: (equation (38)) caused by the neglect of a^"^ will be seen to be of the 

order of 'OOOSttp -: cos 4-6, 

Since the maximimi value of r^jiji^ on this curve is Bc^^, it appears to be legitimate 
to neglect this error. 




Fig. 1. 

Points of Bifurcation on Elliptic Series. 

§ 19. Let us now return to the series of elliptic cylinders, and search for points, of 
bifurcation on, this series. We have found in § 13 that when a^ is finite, but when 
the squares and products of the remaining a/s may be neglected, we have (equation 
(44)) for values of n difierent from. 2, 

fn (a;,, a.2, cig, . . .) = Y — '-^-^-^ -^-''' ^ terms linear in a,^^g, a.^^^^,, &c. . . (61). 

Let us now examine the form, assumed by equation (54) at points upon this series. 
We have, in the first p.lace, df/dcio === 0. 

-ivX /J 
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We have also df^ld/n = whenever m m greater than n. Hence we see that the 
determinant on the left hand of (54) reduces to the products of the terms in its 
leading diagonal^ and that the equation itself is equivalent to the separate equations 



(I 



(X) 



2 



ZjTT p) , 



(62) 



taken for all values of n from 1 to oo , 

Corresponding to a root of (62) there is a point of bifurcation, and the linear series 
starting from the point must be found from, equations (52). From these equations it 
appears that the linear series corresponding to a root of (62) will be such that, as far 
as the first order of small quantities, Og exists only for the values s =. 2, n, n -- 2, 
'}% ~~~" •*,,,, 

Of these series the series n'=^ 2 may be rejected, as corresponding only to a step 
along the series of elliptic cylinders, and not to a new series at all, and the series 
n =: 1 may be rejected, as corresponding merely to a change of origin. 

We are left with the values n ^^ 3, 4, 5, ... , and for any one of these values we 

have, from equation (61), 

3^^ . 

dcL.i, n 1 -~~ 2aa^ 

The points of bifurcation are accordingly given by the equations 

1 1 -f a« 

% 1 ■— 2aa^ 

where /i has the values 3, 4, 5, » . . , 

These points of bifurcation are points on the series of elliptic cjdinders, hence or, 
a.2, and a are connected by equations (57) and (58). If we eliminate w^ and % from 
the three equations (57), (58), and (63), we find, as the equation giving points of 
bifurcation of order n, 



1 1 i- 



m' 



9 



9, 



V, 



(64). 



This equation . must be solved by graphical methods. In fig. 2 the curve which is 
concave to the axis of a is the parabola 




Fig. 2. 
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The remaining curves are the graphs of 

lA ■"""* • « • • • • • « • • • \/ 

for the values n = 3, 4, 5, . . . 

The curve (66) cuts a = at the point y = \\n, and a = 1 at y = 2/n. The 
value of dy/da is a""^. It is therefore obvious that the curves are convex to the 
axis of 06, and since for any value of a the value of dy/da is greatest for that curve 
for which n is least, it is obvious that the curves can never intersect. 

We therefore see that the parabola (65) will meet each of the curves [66) once, 
and once only, for values of a between and 1. Moreover the smaller n is, the 
smaller the value of a at the intersection. 

As we move along the series of elhptic cylinders, the value of a increases from 
to 1. Hence there will be an infinite number of points of bifurcation on this 
series, of orders 3, 4, 5, . . . . The point at which we arrive first is that of order 
n ~ 3 ; those of orders 4, 5, . . . follow in succession. As before^ we find that the 
configuration at the end of the series of elliptic cylinders (a = 1, an infinitely long 
and thin ellipse) is unstable for every vibration. 

The linear series which we expect to be stable is that of order n = 3. To 
find the point of bifurcation of this series we require to solve the equation 
■|- (1 •— c£^) = -g- (1 + ^^) and the solution is found by inspection to be ex = J. 

From equations (52) and (53) we find that at this point of bifurcation co^ = fi^P, 

and ^2 ^^ 5* '^^^ elliptic cylinder at the point of bifurcation is therefore the 

cylinder 

^r, = cr + i{e + 7f) . (67) 

or, in Cartesian co-ordinates, 

g^9. _^ 9^/" = 5a^ .......... (68). 

If we reduce the linear scale of this until the area is a^, we find for its equation 

x^ + 9y^ = 3a^, 

and for its angular momentum, 1*46 times the greatest angular momentum for which 
the circular form is stable. 



Poincare's Series of Pear-shaped Curves, 

§ 20. The configuration of the new linear series of order n = 3 is, near the point 
of bifurcation, of the form 

f>? = ^^ + I {^ + ^^) + h (i^ + V^) + W{^+ V) ' • • • (6^0' 
This new series is seen to be the series corresponding to PoiNCARi's series of pear- 
shaped figures."^" Instead of making a separate problem out of the determination 
of the constants hj and 63, we shall, in order to avoid repetition at a later stage, pass 

^ H. PoiNCARE, ' Acta iMcatli.,' vol. 7, p. 347, 
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at once to the equations determining the general configuration of this series. We 
therefore replace equation (69) by 

§> = a^ + I- (P + 7?2) + 2 &^ ,Co -f 2 ,,C„ (f" + i') ... {70). 

This will be assumed to be the general form of the surface in the linear series 
now under discussion, the quantity 6 being a parameter which vanishes at the point 
of bifurcation. 

The equation expressing explicitly the solution of (70) may be supposed to be 



i). \ 



^^=:(l-,"^-)(-|fo+^\^+f#+f#+...) .... (71), 

\ AlTpJ 

in which £q is the value of f when = 0, and therefore satisfies 

^^-a^ + K^' + r) ..... ^ . . (72)^ 

and f^ is a series of ascending and descending powers of y, say 

^ =r: ao + ajTj + aojf + , . . + ^^--i'^""^ + %^"^ + . . . . . . (73). 

If we calculate the value otYi from (71), we find 

y,=.^^p^n + i^p{i.~~£^TJ,-i-np(i^£^te^u. . . (74), 

where, if ^g is given by (73)^ the value of JJs is 

IJ, = C + a, {i + v) + ia,. (f Jrv')+ ■■■ ■ ■ ■ ■ ■ i^b). 
Using the value of V,- given by (74), the equation to be satisfied at the surface is 

~ ^7? {-rrp - W) + i^P (l - £p) Uo + ^P (l - 2^^) S^"U, := const., 

or, dividing throughout by irp — -Jw^, 

^7) == l^G^JJs + terms independent of ^ ". . , , . . (7C}), 

Equation (76) must be identical witli (70), the right-hand members of both behig 
spherical harmonics, and hence we must have 



H=:l 



and therefore, by equation (75), ^C,, =: cin^Jn for all positive values of n. 

Instead of being given by equation (73), the value of ^^ may now be supposed to be 

given by 

f. ^ .0^ + 2,Go^r] + 3,0^^^ + . . , + a_^frj'"'^- + a^.^^^ + . . . . . (77). 

If we introduce the limitation that the curve is to remain of constant area, we must 
put a„i :== 0. If we now replace a^^, a._3 , . . by new unknowns /J,_]. /-Lg, . . . ,, we 
can write equation (77) in the symmetrical form 



n~ — CO 



in which we know that /L.^ must ultimately be equal to zero, in order that the centre 
of gravity may coincide with the origin (cf, equation 22), 
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Now (sr" will be a function of d, and the relation between a>^ and m at oiir disposal, 
this relation being virtually the definition of 6, We have, however, already assumed 
that at the point of bifurcation ^ ~ and dco^/dd =: 0. We therefore take as the 
relation between co^^ and 6, 

1 - a)y27rp ^ §0 + M^ + §3^'^ + • • • • (79), 

in which Sq ^^ f , and Sg, §3, . . . are as yet undetermined. 
The value of f given by equation (71) is now 

^- fo + 08,e, + 0' (Sof, 4- P.fo) + 0' (80^3 + 8,^1 + fSsfo) + • ■ • • (80). 

If we substitute this value in equation (70), of which we are supposing it to be 
a solution, we obtain 

- C.-3 + 1 1,^3 + ^^, ^ 2i, (^§0^1 + G' (Sof. + fS,^o) + • • •) 

+ t ^^GC,+ s A(e + v'))- • ■ (81), 



A— 1 /; = ! 



and if we equate the coefficients of successive powers of 6 in this, we obtain 

fo^ - «' + I (^0' + ^') (82). 



CO 



UAv~Uo) = fio + tAA^o" + V'1 (83), 



00 



(Sof, + f 8,Q {v ~ Uo) - f (Kiif + 80^1 S n fi,, ^r' 



CO 



+ „Co + t A (fV + t) . . (84)^ 

1 

en CO 



-i 



CO CO 



1 1 

Equation (82) is, as it ought to be, identical with (72). Equation (83) enables us 
to determine the constants which occur in f^. These having been found, equation (84) 
enables us to determine ^o? ^^^^ ^^ on in succession. 

§ 21. As far as first powers of 9, we knoAv that the configuration is of the form 
given by equation (69). We therefore assume at once for ^^ the form 

fi = 3C3772 + Ci - ^ - ^f --• (86), 

in which c,^ is temporarily written for ^C,,. 

Since 8q = f, and since the c's higher then Cg must vanish, we find that equation (83) 
takes the form 

I (>? - 14) ^1 = «3 (fo^' + >?') + <?i (fo + ^) ..... (87). 
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^ , ^^ I ^^^ , 1-000 25,000 

■^'^ A(rf' o^nf o'rj' 



Since ^^ satisfies equation (82), we have 

or, expanding in the appropriate form, and taking a = 1, 

__„ 1 .0.0 

»0 — 1'^ 'T- \ 

From equations (82) and (89) we have 

> O r; > K t '^ I ^1. 

bO -— 2b 07 7 2 -^ 47 I 3 

By a similar process we obtain 

M) — \ 4 ^ 2/ So 2^ 4 7 

Ecf nation (87) can now be put into the form 



1 9 



(88), 



a » 



(89). 



125 



27^3 



2500 

'>r> 4. I * • 



» f 



(90). 



1 •' 

8 7 



■^^ . U X O I o 

4 / -r ^^^ -r g^^:^ 



e » « 



(91). 



8^ 



617 



25 500 

27'??'^ 24-'Sif 



« « $ 






C' 



(? 






00 



V 



4 



» « « 



r , ^ 50 1375 
^3 8^^ + 1^ + a' + ■q^i n + 



■~i 



"^ (]] 






Equating the coefficients of the various powers of 7] we obtain 



1 1/ ^j i« 



8% 



9 



"(/«3 



2^3 1^ 8^''l — 4^3 ! 2^'l 



-2 A/-* 
9 ^3 



6^-1 



8i ^3 









„„ 5 



^ 8 1 "^3 "T 2T^1' 



The first equation is, as it ought to be, an identity. We may assign to Cg any 
value, and therefore take c^ = 1, this being equivalent to fixing the linear scale of 
measurement of #. Solving the remaining equations in succession, we obtain the 
following scheme of values :-— 



a 



Jl< 4 



Gi 



3 > 



C^.i 



0. 



0_^Q 



1 



The vanishing of c, shows that the centre of gravity of the curve is, as it ought to 
be, at the origin. We now have a. the vah,e of f, equation (80), 









3 i^ .>t-:^ j, 1 • • • 



1 



. \ t7 O }» 



§ 22. We now proceed to the determination of 4^3. Equation (84) takes the foiin 



n~oo 



t t « \ / * 



The value of f^ is of the form (equation (78)), 

4- CO 



■00 
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Suppose tills value substituted for $^ in equation (94), and in the equation so 
obtained equate the coefficients of the various powers of '>], For any power of 17 
greater than the fourth it will be seen that equation (94) may be replaced by 



n^' + X! ';?,= » 



t 

The equation obtained by equating to zero the coefficient of any power of rj greater 
than the fourth will therefore be of the form 

a linear homogeneous function of oC„, b^^^^-h? . . . = 0, 

and there is an equation of this form for every value of n greater than 4. This 
system, of equations can only be satisfied by taking 

q\_/„ «»™. oV//-' oV/.'7 — — , , , — — U, 

<v u (■J V (^ I 

We may therefore assume for ^.^ (equation (95)) an expansion of the form 

4 == id^^ + zd^ + 2% + c^i — -j — -^ ;;:' .... (96), 

in which cL is written for oC.. 
From equation (82) we have 

and from equation. (87), 

(^"-|fo)^i-l(^/-l-f)-¥(4 + ^)- ..... (97). 
From these last tAvo equations we obtain 

Ur (3t7 - ¥) - li &H ~ W) ~ ¥^ {Q + y)') + ^F^ (fo + >?)• ■ (98) 

With the help of equations (97) and (98) we can write equation (94) in the form 

Kl oho) V^oS*o -f- oOog-oJ — . 3 2$^! i"- ^1 \ 64^ 487 

-^ ^iv {^^ + ^^) + -^1-^^ do + ^) 

+ <h (^0* + ^'*) + ^4 (fo^'^ + ^^) + ^^2- (4' + ^^J + ^^1 do + v) + 4 • (99). 

It is clear upon examination of this equation that the equations found upon 
equating the coefficients of r/\ tj, t)'^^, &c., will contain only terms multiplied by 

do, c/q, cLx, &c., without constant terms. We therefore take 

.; ^ ./ ^ J _^ _„ A 

Equation (99) now contains only even powers of ^fj. Before we can calculate the 
coefficients of these powers we must obtain series for ||^ and f^/. By squaring 

equation (93) we get 

2000 

9?} 
Next we have from equation (82) 

io^ + Tj^' — -^ {irj — 1). 

Squaring this, and subtracting 2fo^'^^ from each side, 
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^7=z9r,*~-.20^2 + if« + ^^r+' •• ..... (100). 
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^ 4 L ^4 __ 2„5 _ 2 5;^ ^ j lit 3^2 
6o "T" 7 °~^ 4- 2 €^07 t 4- bO 7 » 



Using the series which have already been obtained for ^^ and £q (equations (89) 
and (90)), we obtain 



^0 —16^ T" 6^ ^T 5 4 1 9J0 3 



e can now evakiate that part of the right-hand side of equation (99) which does 
lot contain d,^^ d^, or d^. We have 

S2bl '™- 3^7 8 7 n 7 2 1 i Q_3 H^ 



e / 2 5 5 ^2 3 2 5 \ ._,„.. 7 6 5 ^4 1 1 2 5 ^2 i 1 6 2 5 f ^-^J^^^. I 



e a 3 



34,375 






54r/ 

By addition the sum. of the terms in question is found to be 



(to / O 



45^4 25^2 Lii-7A \. ZZJ^^:^- A^ 

I 6 7 16 7 i 3 6 I g ,^ o [ . . • 

Lastly 3 we have 

Collectin^y thb various series, we find as the form, assumed by equation (99), 



5 25 n 






40/0 



0, 



— 16^ 16^ "T- 3 6 T" g^^ 

Equating the coefficients of the various powers of tj, we obtaui 

4^^'l """ 1 6 t 16^-^4 ? 



- ¥^^4, + f t^3 + m, - - ft + ¥, + 1<^3 , 

2 7 5 ^7 8 0^ B ^7 -^437515000./! 1. i 5_,./ 

Solvingj we obtain in succession, 
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d 



4i 







4.5 
7 ) 

63 



8 



§2. 



QV^g 



d 

d. 



a 



5 6 "T 2 5"%- 



' ^1 



43^75. 

5 4 



992 



81 %• 



We therefore find as the value of l^^ 



*> Hr :"** 



^3 — 7 7 



2. 8. 2 5 ^ 
28 ^ 



4o / i 



0*7 H 



+ ...+S, 



1 984 



(101). 



§ 23. We now proceed to the determination of l^g. Equation (85) takes the form 

+ {S,i,f 3fo + (So4 + mo) (S^o' - ¥) 



CO 



+ 3C0 + t 3C, (^,; + r^") 



(102). 



All the terms on the right-hand side, except those in. the last line, are of odd degree 
in Tj and of degree 5 at most. The same is true of the terms on the left-hand which 
are multiphed by 8^. The terms multipHed by Sg are of even degree^ two at most. 
It is therefore clear that we may at once take 



3^6 



8^7 



3^8 — 3^9 



s'^-^o 



e a « 



= 0, 
= 0, 



2 



e « » e 



V 



ilvo), 



'3 — ; 3O4, — gO^ " 
and assume for ^^ an expansion of the form 

We now calculate the various series which occur in equation (102). We have 

5 25 ^ JOO 

6r} 27??'^ 



So iv -^ iio) ill 



SLnn 

87 



Z4:07}' 



* " ' \ ^' 



3? 



and from equation (97)- 



hiv 



4 
5^ 



)ii^Mi{io' + v')-''f{i. + v): 



64_ 
37597 



8 



N {^<^' + v') + H^'v {So + vn 



This last bracket can be at once calculated from the series of the last page ; we have 

Mv ~ Uo) ii = 8, ih' - (^v + . v^'-\- ■ ■ ■) ■ ■ ■ . (104). 

Next 

fSof, (So4 + i8,Q + (80I, + i8,Q (3^^ ™ i^) =. (fl^ + |S,f„) {U, + 31,3 - J^) 
From equations (93) and (90) v/e have 






125 . 2500\ 



125 . 4000 






. (105). 



From equations (101) and (89) we liave 



N 2 



J ^ 
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8 ^ + 5 ^2M) 



ii 



1 8 0_:3 2 82 5 



8 1 -^T--^ 



s"-^ 






1 q Q, 1 



1 



1 



J^ b 



8 



!'>? + o" + 



80 



)?/ 27'^;'^ 



„-™ 2 2.5_3 



1 4 7 



1_4,]. 2 5 
2 24 / 



21,87^ 






2167?'' 



'|- Uo 



. ^>„ 4. ^«- + l^SO _,. 1 



and hence, by multiplication with (105), 



m 

a / 



2 02 5,^5 



5 6 



-1/ 



I 2 7,1 2 5^3 5_ 8 12 5 
8 9() ^ V n 14 ~ 



78J25 






}- S.> i -!)-T)^ M Qyj ^J^ -^ - _|_ ^ ^ ^ I 



We have from equations (89) and (100) 



4- 







o7] liJ]' 240?]'' 






3 



9r;4 — 20i]" + 'I'"' 



2000 
9^-2 ^ ' * • 



and hence^ by multiplication, we lin^ 



/r /r 2 .-_ 9-,,5 J_ K^,B 10 0.. I 



3500 



81 



V 



5 s e s 



6 2 5 



21,875 



Therefore {Uif i^Q - f 11^^ + Wr - W^ + .l^zr + • 
From equations (91) and (89') "\¥e have 



1 SO^r 3 _ 2825.^ _^ 4 5^,3 1025^^ ZS~I^ l^l^Zll-' JL 

7 ^ bO ""2 8 bo "^ 147 «-- - 5"(y-?J - ^^^^ -|- ]^l^,-'^''^ 



2125 154^75 
84^" '^^ "lUrf 



9 » 



and from equation (93). 



U -ij 'J 






By multiplication of these last two serie.. 



5^ i: /i8 /r 2 2_8.2l/r \ 

^Obl V y bO "^ 2 8 bO/ 



fill 5 
1 1 2 7 



2 6 2.5^3 
"64 V 



li :i 2 5 

3 



r 



50 



We have also 

;^ ^ 4 8 d' /: ™- 6 ^ J 1 ^ U- '' »L "^ JL. 

^0^2 2 5b obi -^ 5 ^2 1 S7 "T o ^ P ,>k. .^ r 

= So {ir + 47? + 0.77- ' + ...} 



«> 8 » 






■¥ + o.7j^3+.,.} 



3 6 7 J ^ "-p . . . 



The last series we require is l'^^. By multiplication of the series (91) and (92), 
we fin 



5000 



t 5 — 15 V 24 3 I 8 7 5 I Z.Z1Z. 4^ 
bo — 3 27 I 48 7 r 2 167 i 9,9^^ t - - 

If w^e now^ collect the various series which have been obtained and substitute them 
in equation (102), we find, as the equivalent of this equation, 
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5 






. _ . __^ 500 

«^ "*" 6>? 27^"' 24V 



e 



^e-;rf + 8637?^^ + e^ ■ ^^ 



+ 83 [1-17'^ — Qrj -\- Q .r) ' + • • 

yY «J-i> rj'T) -'■{— 



Oo 



9. '» 

5 



97; 



2 02 5_5 



.J_ 2 2 5 _ 
l '5 6" 7 



12 7,125^3 f BiJJ 
i\a 71 14 7 



+ 8^ [fij^ -f- 4>^ 4- . >2 '^ + . . . 

78,125 



89 6 



Q 'T K Q 






r • • s 



i 12 8 7 i ( J 4 7 



6 7 5.^5 



•+- Ti 2^ 



2 62 5^3 
" () 4 7 



01 A„ + 21,875 

6250 

V + t; r ■ ■ • 
It] 



3 3 6" 



/ nooo' 



Vf 



+ Co l')?^ + 1^ + q^ + ' • • 

\ 1/ J 



50 
1^ 



5 






7; 



Equating the coefficients of the various powers of t], we obtain 



8 ^5 



6 ^5 I 8 ^'3 



2 r> 



re,+ 



2 500^ 
2 4 3 ^'5 



2 7 ^6 
9 ^3 



I ^3 + 



5.^ 
6^1 



hi ^ 
il20 5^ 

g - Q 



- fii>7 5 I 3 3., 

_ -j_ 2 8" T" S'2^( 

2 ^2 

«~ 16 o 



158.625 125 
896 ^^ 



4 8-^5 "T 8<^3 



'3 



5 6 2A 

28 



8 7r 



|€. + 4^3 + 



3 



r 2 16^5 1" 4^3 ~r 2^1 



.■y 






4 4 0,6 2 5 I 5_0.0 0.^ l^ 5 0^ i 5^ 

"^'"rolfs"" ^ 2 43 % ^^ 9 ^3 i^ 3*^1 



Solving the first two of these equations, we obtain 



^5 



2 2 5 
4 



(110), 



s 



8 6 2 5 

"4 4 g- 



Equations (108) and (109) may be written 



6 2 5^ 
7 2 ^5 



J 5 

4 



'i^B + To^i)- 168^ 



56 2 5 
28 



21 00 ^ 

~ Hi H 



2 1 
3 



^ {en + A^i) 



6 2 0^ 

9 ^^ 



8^1 



4 4 0,6 2 J 
To 8 



(lOG), 

(107), 
(108), 
(109). 



(111). 

(112), 
(113). 



Multiply (112) and (113) by 20 and 9 respectively, and subtract, and we obtain 

- 03 7_5 h-\A\ 

From eqxiations (110) and (ill) we obtain 



H^e, + 3008, 



X^p 

2 ^-^-1 



6 2 5 



9 3 7 5 



and hence, by comparison with (114), e^^ =: 0. 

This vanishing of e_i suppHes a searching test of the accuracy of the work. 
Equation (112) becomes, after simpHfication, 

^3 n^ 10^1"-" 168 • • • • • • • • VLij;. 

Our equations do not enable us to determine e^ and ej^ separately ; they are, 
however, all satisfied by taking 

e^ r= "— ~i^ X « . . ( 1 J. / ) 



e 



17 7 5 



16 8^ + X • • • (11^)5 
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m. 



where X is unknown. We therefore have (equation 103) 

^3 ==^ lip,/' ^ i^^3^^-77^ + X (37?^- ^^^- ^^f) + terms in ri^', ^q~\ &c. . (118). 

Using the valne of So. given by equation (111), we find as the values of d^, d^^ 
and d^ (p. 91), 



d 



4-5. 
7 ? 



9 5 5 



d 







6 9 8 7^ 5 
"50 4 



§ 24. Collecting the vahies of the various constants, we find as the equation to 
the surface (equation 70), 



+ ^' m^ ie + v') - - flF (f '^ -!- V"")} + terms in e\ 6% &c. . (119). 
The occurrence of the indeterminate quantity X can. easily be accounted for. For 

if we have a solution 

f>? = «2 + Of, + 0% + ^!/;, + . . (120), 

corresponding to a parameter 6 which is connected with the rotati.on by the relation 

1 -^ (^1217 p = So + S,.^-^ + S.ff^ + ... , . , . . (121), 

then we can obtain this samB solution in another form by replacing the parameter 6 
by a new parameter 6 + X^^. As far as 6^ this leaves the relation (121) betv/een cu^ 
and 6 unaltered^ whilst the equation to the surfaces as far as ff^ becomes 

^'Yj^fl-i^ie + W)j\-V6%r\r 0%+ . . . „ . . . (122). 

It accordingly appears that in equation (119) the value of X is entirely at our 
disposal We shall therefore take X ^^ 



'nvestigation of tkability. 

§ 25. There is a large d priori probability that the linear series w^e are now 
considering wdll be stable for small values of 9, but it will be well to rigorously 
examine the question. 

It appears from PoiN care's researches that the whole investigation reduces to 
determining whether the angular momentum is a maximum or a minimum at ^ = 0."^ 
We therefore require to calculate the angular momentum as far as 6^, and the answer 
to our problem wdll depend upon the sign of the term containing 6^, 

As far as 6^^ the equation to the surface in polar co-ordinates (i\ <^) is 
(equation (119)) 

r- = 1 + |-r'- cos 2(^ + 2Q {f cos 3^j.— ^f-r cos ^y) 

^ (^r^ cos 4^ ^- ^^1^ cos 2c^ + ^^Wr) ^ ^ ^ ^ . (123). 



The moment of inertia is given by 

I 



r^' dv r d<j) 



"^ H. P01NCARE5 "Sur la Stabilite de FEquilibrc des Figures Pyriforines 
vol. 198, p. 333. 



" I T) 



9 8 » ^ 



X llll. xTRllS.^ l\% 
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taken over the surface^ arid therefore by 

.L ^^^^^ "o" "^ ' Cv(p • . . . . . . . . . il^-tcK 

J 

where r is given by equation (123). 

Let us assume a solution as far as 0'^ of the form 

r'^=za^l + /3e + y0^) .,.,.... (125), 
where a, /3, y are functions of c^ only, a being given by 

fx'^(l — f COS 2(f)) =:!...,...,. (126). 



The value of ?*^ corresponding to solution (12 5) is 

and if we substitute this value for r'\ and the similar vahies for r and r^ in equation 
(123), and equate the coefficients of 9 and 6''^, we obtain 

-|-^ =: 2 (a'^ cos 3cp — -£-oc cos c/>) ../,,, (127), 

iy _ ps ^ fcc^^ COS 3(^ --- |a^ cos ^ + -^a^ cos 4(^ »- ^f -^a^ cos 2(^ + ^f f P, 

whence, by elimination of /3, 

|-y zz: 8 (<x'^ cos 3^ -^ i^a cos (/>) (a^ COS 3<j?) — fa cos 6) 

-j- ~->^Ci COS 4(p -^' ' q — Ct COS -djcp -f- — g"Q-^"~. 

We can eliminate c/) from tliis equation by the help of equation (126). The 
resulting value for |-y contains only even powers of a ; if we simplify this, and 
transform the numerical coefficients to decimals, we tind 



iy — m-la}' — 319^8^* + 264'4a^ + 159*4. 

ow we require to find the coefficient of ^'^ in I (equation (124)), and this is 

M a^y dip. We therefore require to know the value of a^^^d(j} for n = 2, 3, 4, 5. 

Jo Jo 

This integral can easily be evaluated for all positive integral values of n ; the values 
which we require at present are as follows :— 



1 f^ 
TT Jo 



4-6 17-0 70-1 305-3 

TT 



■TT 



Substituting these values, we find at once that ^ a^y dcf) is a positive quantity. 

J 
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Thus the moment of inertia and the angular velocity both increase as 6^ increases from 
the value ^ = 0. The moment of momentum is therefore a minimum for the value 
^ z=: 0, and this proves the stability of the series of pear-shaped figures. 



The Seines of Pear-Shaped Oarves, 

§ 26. The equation to the curves of this linear series has already been calculated as 
far as 6^, In order to obtain a still better idea of the shape of the curves I have carried 
the calculation two degrees further. The calculation of these last two degrees is 
extremely heavy, and I have omitted all details in order to save space. The method 
is precisely similar to that which was followed in the calculations of §§20-23. 

It was found that the coefficients multiplying terms in 6^ and ^^ w^ere inconveniently 
large, and to obviate this, the parameter has been changed from ^to 10^'^^. After 
making this change we find, as far as ^^, for the equation to the surface expressed in 
polar co-ordinates, and for the equation determining w^, 

r^ = (1 + -139^^ + -02319^ + _ .) — =211^r cos(/» 

+ (-8 - 138^'^ — -069^^ + . _) r^ cos 2^ 

+ ('0636^ — '0064<9''^ - '00316^^ , . .) r^ cos ^ 

+ (-0136^^ + '0008^^' + . . .) r"" cos 4^ 

+ (-0036^^ + •00093^9^^ . . .) r^ cos o^ 

+ (•0011^^^^+ , . ,) r«cos (k/> 

+ ('00043^^ + . . .) T^ cos 7(j^ + &c. . . . (128), 

1 — or/27Tp = '625 - '0196' — '01619^' (129). 

§ 27. We must next consider within what limits the calculated terms of equation 
(128) will give a good approximation to the complete equation. It is clear that for 
given values of r and 6 the worst approximation may be expected when <^ :=: 0. 
Let us therefore consider the function ^ (r, 6)^ defined by 

* (r^ e) ~ (1 + 139^^^ + -0236^'' + . , . ) - ^211 ^r 

- (-2 + 138^^ + •069^^' + . , . ) 1- 

+ (-0636^ — =00646>3 »»- •OO31^^0r^ 

+ (*01 3^2+ '0008^^^+ . . . )r'^ + ('0036^^+/00093^^)r^ 

+ ('0011^'^ + . « • ) r^' + ('00043^5 ^ _ ) ^.7 + &c. . . (130). 

The value of <E> (r, 6) is expressed by a doubly infinite series, of which only a few 
terms are known. When r = 0, ^ == 0, the value of <J> is known to be accurately 
equal to unity. For small values of r and 6, equation (130) will give <I> with 
considerable accuracy, but for larger values of r and 6, the terms calculated will be 
inadec|uate to give a good approximation to the value of #, What then, we 
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inquire, are the values of r and 6 over which this approximation may be regarded 
as good ? 

The coefficient of each power of r is an infinite series of powers of 6, of which 
all terms up to 6^ have been calculated. A glance at these series will show that the 
approximation is tolerably good so long as 6^ <. I, but begins to break down as soon 
as 6 exceeds this unit value. 

Supposing that we have assigned to d some definite value less than unity, the 
value of *(r, 6) will be given by an infinite series of powers of r, of which only 
the first seven are known. For small values of r these first few terms will give a 
sufficiently good approximation, for larger values the approximation will be bad. 




Fig. 3. 



and for still greater values the series will become divergent^ so that the first few 
terms give no approximation at all. It will be seen from inspection of equation (128) 
that the approximation will be tolerably good so long as r^ < 1/^^. 

The conditions under which the calculated terms will give a good approximation 
may accordingly be supposed to be that ^^ < 1, and r^ < IjO^, In fig. 3 is 
represented the plane of r, 6, The part of this plane over which the approximation 
is good is that boimded by the four curves 

This is the portion which is shaded in the figure. 
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In this same figure the thick curves represent the locus, 

# {r, 6) = 0, 

calculated upon the supposition that the calculated terms oJ: $(?% 6) give a 
sufficiently good approximation to the whole. For the greater part of the curve 
this assumption is not justifiable, so that the curve requires adjustment, the 
amount of this adjustment increasing as we recede from the shaded portions of the 
plane. The most important points on the curve are those at which dd/dr = 0. 
These may with sufficient accuracy for our present purpose be taken to be r = 2, 
=2 l^ and r =:: — 2, ^ =:: — » 1 . 

§ 28. The points at which the axis (^ = meets the curve of which the equation is 
(128) are given by 

<|) {r, 0) — 0. 

Fig. 3 accordingly enables us to trace the motion of these points as we move aloug 
the linear series, i.e.^ as increases from, zero upwards. At ^ = we have, of course, 
two equal and opposite roots — r ™ i y^5. As increases the positive root increases, 




JL A Wi * ^- e 



while the negative root numerically decreases. Remembering that the centre of 
gravity of the curve must remain at the origin, we see that this indicates a general 
thickening of the half of the curve in which ^ > 7t/2, with a diminution in the 
thickness of the forward half, and consequent lengthening of this half These 
features become more marked as increases, until we reach the value :=^ 1, at which 
a new feature presents itself For here there are two new roots occurring at the 
point r = 2. This indicates that the fluid separates into two portions when the 
value ^ =: 1 is reached, the point of separation being r = 2 (approximately). We 
are at once struck by the great inequality in size between the primary and satellite : 
the former extending approximately from ?• = — 2 to r = + 2, and the latter only 
from r = 2 to r = 2|-, The ratio of the linear dimensions will therefore be some- 
thing like 8 to 1 , but it must be remembered that our results require considerable 
correction on account of the imperfections in our approximations. 

§ 29. In fig. 4 the thick curve is the elliptic cylinder corresponding to ^ = 0, and 
the dotted curve is the adjacent curve corresponding to a small value oi {0 ^ |). 
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In figs. 5 and 6 the curves are those corresponding to 6^ = -^ and ^^ = 1. A glance 
at fig. 3 will show that there are difficulties in the way of drawing these latter curves 
with much accuracy. I have given in detail some of the calculations used in 
drawing the curve 6^ = 1, in order that the reader nmy judge for himself as to the 
closeness or otherwise of the approximations. The curve ^^ = l is of course much 
easier. Before passing on to the calculations, two points ought to be noticed. 

(i.) It will be noticed that in the various ^-series (the coefficients of powers of r in 
equation (128)) the terms last calculated are without exception of the same sign as 
those previously calculated. There is therefore some justification for hoping that 
the remainders in these series will be of the same sign as these last terms. If this is 
so, the error introduced by the neglect of these remainders could, to an appreciable 
extent, be reduced by an adjustment in the value of 0, Thus we shall be attempting 
to calculate the curve for (say) ^ = 1, and shall obtain a curve which is much more 
like the curve for some smaller value of (say =■ '98) than it is like the curve 
^ = R Regarded as an attempt at tracing a surface of equilibrium the error will be 
much smaller than if regarded as an attempt at tracing the particular curve =. 1, 




Fig. 5. 

(ii.) It will be noticed that the sign of the leading terms in each of the series 
multiplying r*, r^, r®, r'^, is positive. An examination of the method by which these 
leading terms are calculated will show that this is a general law : all the leading 
terms after r^ are of positive sign. Thus the error will be reduced by supposing the 
series (128) continued to higher powers of r by a suitably chosen series of terms. 
I have accordingly done this in the calculations, and the conjectural terms are, 
throughout, enclosed in square brackets/^' 



The Curve 0^ 



2 



(Figo 5.) 



§ 30. In tracing this curve I started from (j) = rr, and calculated a series of points 
on the curve for decreasing values of ^. The approximation at ^ = tt was found to 



•X- 



The efi'ect of these corrections must, of course, be small ; but, at the same time, it seems as well to 
make use of any definite knowledge that we possess, 

o 2 
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be good, the root being 1*98, and the error occurring only in the third decimal place. 

As cj) decreases the approximation improves, and at </> = 7r/2 the error occurs only 

in the fifth decimal place. At ^ ~ 7|-^ the error again appears in the third place, 

and after this the approximation is bad. My plan was to calculate for smaller values 

of (f) as well as I could, taking care to keep the values of r in defect rather than 

excess of their true values. The curve was then plotted out on i)aper ruled with 

squares of 1 millim., the unit of length being taken to be 50 millinis.''^' 

The area of the elliptic cylinder of fig. 4 is known to be 

13,090 sq. millims,, 

and this would also have been the area of the present curve had it been accurately 

drawn. The area of the curve (obtained by counting squares) was, however, found 

to be 

12,776 sq. millims. 

The moments about the axis (]j = 7r/2 of the two parts of the curve {(j) > 7r/2 and 




Fig. 6. 

(/)<7r/2) ought of course to be equal : these w^ere found to be respectively 

298,290 cub. millims. and 302,850 cub. millims. 
It was therefore obvious that the curve had been too much shortened in the region 

in which ^ < 7|^°. 

Readjusting the curve in this region so as to divide the error as equally as possible, 
I arrived at a curve with outstanding errors in area and moment of 

139 sq. millims. and — 139 sq. millims, at r := 125 millims., ^ = 0. 

This is the curve given in fig. 5. It will be seen that the error is one of about 
1 per cent. 

Calculation of the Curve d =-- 1. (Fig. 6.) 

§ 31. The equation of the curve is found to be 

r^ — 1-162 — -21 Ir cos(/> + '593?^^^ cos 2(j& + •053r^cos 3</) + 'Ol^r^ cos 4()f> 
+ •0045r' cos b(j> + -001 Ir^ cos ^ + •00043r7 cos 7^ 
+ [-0002^^008 8^ 4- •0001r^cos9^+ . , .] 
the terms in square brackets being those mentioned at the end of § 29. 

■^^ This has been photographically reduced to 25 millims. before printing. 
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W hen (j) = n, * 



When r = rS , <[> 



fM^/dr 



therefore 
When <f> 



■001 Ir" 



•001] 



We calctilate in succession the value of the radius vector when <^= 180°, 150°, 120°, 
90°, 60°, 45°, &c. The value of r is a root of <I>(r) = 0, where 

cp (r) = — r^ + ri62 — ■211r cos<)f> + &c. 
Our method will be to find a value of r for which * (v) is small, and calculate <J? (r) 
and cl^jdr at this point. The true root R is then given by 

(M3 {r)ldr ' 
1-162 -f •211r — ■407'r2- -0537^3 + •014r* — •0045r3 + -OOllrs 
~ ■00043r7 + [-0002^8 - "000 Ir" + . . /. 

= 1-162 + -380 - 1-319 — -310 + -147 - "085 + -038 — '026 
+ [-020 - -018 +...] = -- -003 + [-01], 

: -21 - 1-4 — -5 + -3 — -2 + -1 — . . . = — 1-4, 

R = 1-8 - -002 + [-007], or, R = 1-805. 

57r/6, 't = 1-162 + -182^ - •704r3 - -OOr/^' + -OOSgr^ 

H- -OOOSer" - [-000 Ir^ + -OOOOr^o +...]. 

When r=v''2 = l "414, # = l-162 + -257 — 1-408--0284- -022 — -009 + -004 

= -000 ~ [-001], 
therefore R= 1-414. 

When (^ = ^/2, # = i-j 62 - l-593r^ + -OUr* - -OOllrS + [•000-2r8] + . . . 
When r^ — '74:, r = -859 ; 

'P = 1-162 - 1-179 4- -007 ~ -0003 + [-0000] + ... = - -QIO ; 
d^/dr = — 2-5 ; therefore R = -859 — -004 = -855. 

When <f} = 7|-°, * = 1-162 -- •209'r - -427r-2 + -046r3 4--012r* 4--0036rS + -OOOSr^ 

+ -OOOSr^ + [-OOOlrS + -00004r9 -|- . . .]. 

When r^ = 3, r= 1-732, <1> = 1-162 ~ -348 — 1-281 + -241 + -108 + -056 + -021 

+ -014 + [-008 + -005 + ...]=- -027 + [-01]. 

- 1-44 + -42 + '23 + '17 + '08 + . . . = — '75 ; therefore 
R = 1-732 -- -036 + [-01] = 1-71. 

: 1-162 — -21 Ir — -410r2 + -052r'^ + -0141'* + -0047^ -f -OOlr^ 
+ -0004r'^ -f- [-00017r8 -f -OOOOSr" + . . . 

1-162 ~ -422 — 1-640 4- '416 + "224 + -128 + -064 + -051 
+ [-05 + -04 + ...]=- -02 + [-1 ?]. 

- 1-640 + -624 + -448 + -320 + -192 + -178 + ['18 -f -18 +. . .] 
=. -01+[?]. 

= 1-162 ~ -390 - 1-328 + "297 + -147 + -076 + -034 + -023 
+ [-017 + -014 + ...]= -02 + r-05]. 

- 1-47 + -49 + -33 + -21 4- -11 + -09 + [-08 + -07 + . . .] 
= - -45 + [-2]. 



d^/dr = — -209 



lien f =: 4 , # 



When r = 2, # 



d^/d 



Ir 



•211 



When r = 1*8, # 



d^P/dr 



•211 



102 



MR. J, H. JEANS ON THE EQUlLIBRltJM 



There is therefore a root at about r = 1'95, but it is clear that we are ah-eady in 
the region at which the approximation ceases to be satisfactory, and that we are close 
upon the region in which the series become actually divergent. 

The smallness of d^jdv indicates that somewhere in the neighbourhood of the point 
just calcula^ted we come to a point at which there is a pair of equal roots for r, and 
therefore a ^'minimum" in the value of ^. This is the ^^neck" of which the first 
signs are apparent in fig. 5. Let us refer to all the matter to the left of this '^ neck " 
as the '^ primary/' to all that to the right as the " satellite." Let the exact line of 
division be a vertical line at a distance 2 from the origin. 

The primary has been drawn with fair accuracy ; the satellite must be drawn in 
the manner adopted in the difficult region of the former curve. 

The area of the primary was found to be 

11778 sq. milhms.3 

and this leaves 1312 sq. millims. to be accounted for by the satellite and the error in 
drawing. The centre of gravity of the primary was found to be 4*7 millims. to the 
left of the origin. Distributing the error as equally as possible, I have arrived at the 
curve of fig. 6. The area of primary and satellite are respectively 

11778 and 881 sq. millims., 

the error in area is a defect of 431 sq. millims,, and that in the moment about (^~7r/2 
is that of an excess of 431 sq. millims. at the point r = 125 millims., <^ = 0. The 
error in the whole curve is therefore about 3 per cent. ; that in the satellite is 
unfortunately of the same order of magnitude as the satellite itself 

§ 32, The following table sums up the results which have been obtained, and also 
contains some new^ results. The moments of momentum of the last curves were 
obtained by a process of coiuiting on squared paper, and are not carried to any great 
accuracv. 





I 

Curve. 




Area. 


(1)2 

27r^o 


Angular momentum 
( omitting factor A/ ^ 1 . 


(1) 

(2) 
(3) 
(4) 
(5) 

(6) 

(7) 
(8) 

(9) 
(10) 


Circle .... 

5, ... a 

,5 .... 

Ellipse .... 
„ fig. 4 . , 

J^'ig. 5 . . . . 

Fig. 6 . . . . 

,, .... 

1 The two parts 
J of curve (8) \ 


Point of bifurcation , . 
Ratio of axes ^/5 : 1 . . 
Point of bifurcation : ratio 
of axes 3 : 1 

^ = \ » • 

^2:^1 . . . ... . 

Separation of fluid into 
primary and satellite 

Primary 

Satellite 


1 
1 
1 

1 
1 

•931 
•0711 



'43 

•43 

'375 

•39 

•421 
'431 

•431 

•431 



•33 
•35 
•44 
•51 

•53 

•551 
•571 

•401 

•1711 
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In this table the quantities of which the numerical values are doubtful are 
marked with a query. A single query indicates a probable error of 1 or 2 per 
cent. ; a double query indicates that the error may be comparable with the quantity 
itself. 

Let us examine the state of things just after separation has taken place. The 
satellite is describing an orbit about the primary, both bodies rotating with the same 
angular velocity. This angular velocity is, to within a few per cent., given by 

o^V'^Trp =: '43 . ......... (131). 

If the satellite exerted no attraction upon the primary, the figure of the primary 
would be a figure of equilibrium under the influence of a rotation given by (131). 
The force exerted by the satellite may be divided into two parts, a iniiform force 
in the direction :=: 0. and a tide-generating force of the usual kind, if the former 
of these existed alone, the configuration of the primary would still be one of 
equihbrium under a rotation of amount given by (131). We therefore see that the 
actual configuration of the primary may be regarded as a configuration of equilibrium 
under rotation given by (131), disturbed by the tide-generating potential which is 
caused by the satellite. 

Since this tide-generating potential is small, except in the immediate neighbour- 
hood of the satellite, it ought to be possible to remove the tides from the surface 
of the primary, and form a pretty good idea of the configuration which would be 
the configuration of the primary except for tidal disturbance. If this is done with 
the primary of fig. 6, it will be found that the remaining curve is a very good 
ellipse. We may therefore conjecture that curve (9) is an ellipse deformed by the 
tidal influence of its satellite. 

Now the ellipse corresponding to the amount of rotation given by (131) is 
curve 4 of the preceding table. We see that the axes are in the ratio y^5 : 1, and 
this is in good agreement with the ellipse obtained by removing the tides in fig. 6. 
The momentum of the elhpse of unit area of which the axes are in the ratio .y/5 : 1 
(curve 4) is '44. If we reduce this so as to apply to an ellipse of area '93 instead of 
to one of unit area, we find an angular momentum of '38. Since this ellipse must be 
supposed to rotate not about its centre, but about the centre of gravity of itself and 
a satellite about one-fifteenth of its mass, situated at the end of its axis, this angular 
momentum must be increased to about %39. The small discrepancy between this and 
the value '40 obtained for curve 9 may be accounted for partly by errors of 
approximation, and partly by the increase of mouientum caused by the tidal 
deformation of the ellipse. 

We can check our result in another way. The equation of the ellipse being 

ax^ + hy^ = 1 (132), 

the force at ci:^, 0, a point near the extremity of the major axis and outside the ellipse, 
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may, to a good approximation, be taken to be -, — —— -• ."" ^, and if x^ is the 
centre of gravity of the satellite^ this mnst be equal to 

the factor '93 being introduced to allow for the displacement of the centre of gravity 
of the primary. 

Taking y^/h = .y^/{pa)^ we find the equation 



^ ^ '93 " 



n 



and putting oj^'/27rp =^ '43, this gives the values 

aotr :=: 1*54, ^Vc . x ^ 1°24. 

Now ^ya . X is the sum of the semi-axes of primary and satellite divided by the 
semi-axes of the primary. The equation just found is therefore about as true as 
could be expected, the linear diameters of primary and satellite being approximately 
in the ratio of 4 to 1. 

The ellipse given as curve 4 is stable, and, since the mass of the satellite is small 
compared with that of the primary, we may suppose the combination of primary and 
satellite to be stable. Thus, if our conjecture as to the interpretation of curve 8 
is sound, it appears that the linear series which commences with curve 5 remains 
stable until the mass separates into two masses. 

The motion of a gradually-cooling mass will therefore be through the following 
cycle of changes. Firstly, inci'ease of the ratio [angular momentum -— (area)^] until 
we reach curve 5. Then motion along Poincaee's linear series until we reach 
curve 9. At this point separation takes place, and the primary is left as a tidally- 
distorted form of curve 4. As the satellite recedes the tidal distortion decreases, 
and as the value of the ratio [angular momentum -^ (area)^J again increases, the 
configuration moves along the Jacobian series of elliptic cylinders until curve 5 is 
again reached. This completes the cycle, and the continual repetition of the cycle 
can only be ended by solidification, or some similar cause which is outside our present 
considerations. 

I have not attempted to give any discussion of results from the point of view 
of dynamical astronomy. The complications introduced by the heterogeneity and 
compressibility of natural substances, as well as by the difiference between the two- 
dimensional and three-dimensional problems, are so great that any discussion with 
reference to the actual conditions of astronomy would be impossible in the present 
paper. 

I have had the advantage of frequent conversations with Professor Darwin on 
the subject of this paper ; my thanks are also due to Professor Forsyth for advice in 
connection with the earlier sections, 



